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Inclusive transverse momentum spectra of r\ mesons have been measured within pj = 2 - 10 GeV/c at mid- 
rapidity by the PHENIX experiment in Au+Au collisions at ^JsT^ = 200 GeV. In central Au+Au the T) yields 
are significantly suppressed compared to peripheral Au+Au, d+Au and p+p yields scaled by the corresponding 
number of nucleon-nucleon collisions. The magnitude, centrality and pj dependence of the suppression is 
common, within errors, for T) and Jt°. The ratio of T) to it spectra at high pj amounts to 0.40 < R-^/rft < 0.48 
for the three systems in agreement with the world average measured in hadronic and nuclear reactions and, at 
large scaled momentum, in e + e~ collisions. 

PACS numbers: 25.75.Dw 



The major motivation for the study of high energy nucleus- 
nucleus (A+A) collisions is the opportunity to probe strongly 
interacting matter at extremely high energy densities. Of par- 
ticular interest are energy densities well above the expected 
transition from normal hadronic matter to a deconfined sys- 
tem of quarks and gluons. Lattice Quantum Chromodynamics 



(QCD) calculations Qt predict that this transition will occur at 
a temperature of T w 170 MeV 10 12 K. The formation of a 
quark-gluon plasma (QGP) in A+A reactions should manifest 
itself in a variety of experimental signatures \ 2]. 

At center-of-mass energies reached by the Relativistic 
Heavy Ion Collider (RHIC), arguably the most exciting ex- 
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perimental results so far are connected with the predicted "jet 
quenching" phenomenon |3j|4j,|5|] due to energy loss of hard- 
scattered partons as they traverse the dense medium formed 
in the reaction. Since (leading) hadrons with pj > 4 GeV/c 
at RHIC carry a large fraction of the momentum of the parent 
quark or gluon ((z) = Phadmn/ Pparton ~ 0.5 - 0.7 0,0]), par- 
ton energy loss results in a significantly suppressed produc- 
tion of high-/?T_hadrons |4J]. The inclusive spectra of high-/^ 
neutral pions 1 8, 9] and charged hadrons II (ill ill in Au+Au at 
y/s NN = 200 Ge V are indeed suppressed by as much as a factor 
of five compared to the corresponding yields in p+p fl2il and 
d+Au II 3lll4l . scaled by the number of incoherent nucleon- 
nucleon (AW) collisions. The centrality IU5I1 . p T Plnlllill 
and center-of-mass energy 1 19] dependences of the observed 
quenching are consistent with theoretical calculations of QCD 
energy loss due to multiple gluon emission in a dense medium. 
Assuming a thermalized parton system, the magnitude of the 
suppression for central Au+Au at y/s NN = 200 GeV implies 
initial energy densities above 15 GeV/fm 3 , ^100 times larger 
than normal nuclear matter lEoll . 

The equal amount of suppression for 7C° and /z* observed 
above pj w 5 GeV/c for the same Au+Au centrality seems 
to indicate that the mechanism of quenching is independent 
of the identity of the high-/?^ light-quark hadron. This is 
expected if the suppression takes place at the parton level 
prior to its fragmentation into a given hadron. Indeed, in 
this case the high-p^ deficit depends only on the energy lost 
in the medium by the parent (u,d,s) quark or gluon and not 
on the nature of the final leading hadron which will be pro- 
duced with the same universal probabilities (fragmentation 
functions) which govern hadron production in the vacuum 
in more elementary systems. The partons involved in high- 
pr hadroproduction considered in this work have typical mo- 
menta > 5 GeV/c, ten times larger than the "bulk" average 
momenta (pj) ~ 0.55 GeV/c of the system |21]. Such ener- 
getic partons are then supposed to traverse (and lose energy 
in) the medium and hadronize in the vacuum a few tens of 
fm/c later II 1 311 . The equal suppression of ft* and 7C° does not 
by itself provide a conclusive argument for parton energy loss 
before fragmentation in the vacuum because above pj ~ 5 
GeV/c, unidentified charged hadron yields are dominated by 
TC* II ill . Measurement of the yields of an additional light- 
quark species like the r| meson at large enough pj allows a 
confirmation of the independence of the quenching with re- 
spect to the nature of the produced hadron, and tests the con- 
sistency of the data with medium-induced partonic energy 
loss prior to vacuum hadronization. Besides its interest as a 
signal in its own right, the r| meson constitutes, after the 7t°, 
the second most important source of decay e and y. Reliable 
knowledge of their production cross-sections is thus required 
in order to eliminate the background of secondary e and y in 
single electron J22I, dielectron |23] and direct y 1E5I1 measure- 
ments. 

This Letter presents measurements of the r| meson by the 
PHENIX experiment 12611 in Au+Au collisions at y/s NN = 
200 GeV during the second RHIC run (2001-2002) and com- 



pares them t o ri from p+p and d+Au and to 7C° § 1 
and direct y |25] from Au+Au, all measured in the same ex- 
periment at the same ^/s~^. The r| measurement reaches the 
second largest pj for identified hadrons at RHIC, after the 7C°. 
The analysis reported here uses Beam-Beam Counters (BBC, 
3.0 < I Ti| < 3.9) plus the Zero Degree Calorimeters (ZDC) for 
trigger and global event characterization. For each collision, 
the reaction centrality is obtained by cuts in the correlated dis- 
tribution of the charge detected in the BBC and the energy 
measured in the ZDC 12811 . A Glauber Monte Carlo model 
combined with a simulation of BBC and ZDC responses is 
used to determine the corresponding nuclear overlap function 
(Taa) for each centrality |8]. The r| mesons are reconstructed 
at mid-rapidity in the lead-scintillator (PbSc) electromagnetic 
calorimeter 1 29] via their yy decay mode (BR = 39.43%). The 
PbSc consists of 15,552 individual lead-scintillator sandwich 
modules (5.54 cm x 5.54 cm x 37.5 cm, 18 Xo), grouped in 
six sectors located at a radial distance of 5.1 m from the beam 
line, covering a total solid angle of Ar| w 0.7 and A0 w 3tc/4 
rad. The energy calibration of the PbSc modules is obtained 
from the beam-test values and confirmed with the measured 
position of the 7C° mass peak, the energy deposited by mini- 
mum ionizing particles traversing the calorimeter, as well as 
with the expected Epbsd Ptmcking ~ 1 value for e identified by 
the Ring-Imaging Cerenkov detector. The systematic error on 
the absolute energy scale is less than 1.5%, which translates 
into a maximum 8% uncertainty in the final r| yields. 

For this analysis a minimum bias (MB) trigger sample of 
34 x 10 6 events, also used for the previously published 7C° 
analysis |8], is combined with a Level-2 trigger event sample 
equivalent to an additional 30 x 10 6 minimum bias events 
The Level-2 trigger sample is obtained with a software trig- 
ger on highly energetic particles (3.5 GeV threshold). The 
resulting trigger reaches a 50% (100%) efficiency forr| above 
Pt =5(7) GeV/c. The normalization of the Level-2 data 
sample relative to the MB data sample is accurate to 2%. Both 
sets of events are required to have a vertex position |z| < 30 cm 
along the beam axis. Photon candidates are identified in the 
PbSc by applying particle identification (PID) cuts based on 
the time-of-flight and shower profile The systematic 

uncertainty on the yields related to the applied PID cuts is 
~8%. The T| yields are determined by an invariant mass analy- 
sis of photon pairs with asymmetries \Ey\ —E~Ji\/(Ey\ +E^) < 
0.5. The combinatorial background is obtained by combining 
uncorrelated photon pairs from different events with similar 
centrality and vertex, and by normalizing the distribution in a 
region below (m inv = 400 - 450 MeV/c 2 ) and above (m !m , = 
750-1000 MeV/c 2 )theri mass peak. The resulting distribu- 
tion is fit to a Gaussian plus an exponential to account for the 
residual background not described by the mixed event back- 
ground (inset of Fig. The open (solid) symbols depict the T| 
signal after mixed (plus residual) background subtraction. To 
estimate the uncertainty in the subtraction procedure, different 
pair asymmetries and an alternative linear parametrization of 
the residual background are used. The signal-to-background 
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FIG. 2: Nuclear modification factors for T| in Au+Au centralities: 
0-20%, 20-60%, 60-92%. The error bars show point-to-point un- 
certainties. The absolute normalization error bands at Raa = 1 show 
the uncertainties in (Taa) for decreasing centralities. The error box 
on the right indicates the 9.7% p+p cross-section uncertainty fbUl . 



FIG. 1 : Invariant r\ yields as a function of transverse momentum for 3 
centralities and MB Au+Au at y/s NN = 200 GeV scaled by the factors 
indicated in the plot. Inset: Invariant mass distribution of y pairs with 
pr =4-5 GeV/c measured in MB Au+Au, after mixed-event (black 
open circles) plus residual (red solid circles) background subtraction. 



ratio in peripheral (central) collisions is approximately 1.3 
(1.5) for the highest p T and 0.05 (0.002) for the lowest p T . 

The raw spectra are normalized to one unit of rapidity 
and full azimuth. This purely geometrical acceptance fac- 
tor amounts to ^4 at large pj. The spectra are further cor- 
rected for the detector response (energy resolution, dead ar- 
eas), the reconstruction efficiency (analysis cuts), and occu- 
pancy effects (cluster overlaps). These corrections are quan- 
tified by embedding simulated single r| from a full PHENIX 
GEANT simulation into real events, and analyzing the 
merged events with the same analysis cuts used to obtain the 
real yields. The total r| yield efficiency correction is ~3 and 
rises <20% with centrality. The losses are dominated by fidu- 
cial and asymmetry cuts. The nominal energy resolution is 
adjusted in the simulation by adding a p j -independent energy 
smearing of 3% for each PbSc tower. The shape, position, and 
width of the r| peak measured for all p^'s and centralities are 
well reproduced by the embedded data. 

The main sources of systematic errors in the measurement 
are the uncertainties in the yield extraction (10-30%), the 
yield correction (10%), and the energy scale (a maximum of 
8%). The final combined systematic errors on the spectra 
are at the level of 10-30% (point-to-point) and 10-20% (p T - 
correlated) depending on the pj and centrality bin 1E7I1 . A 
correction in the yield to account for the true mean value of 
each pj bin is applied to the steeply falling spectra. The fully 



corrected pj distributions are shown in Fig^for MB and 3 
centrality bins (0-20%, 20-60% and 60-92%) scaled for clar- 
ity by the factors indicated. The error bars are the quadratic 
sum of statistical and systematic errors. 

Medium effects in A+A collisions are quantitatively deter- 
mined using the nuclear modification factor given as the ratio 
of the measured A+A invariant yield over the p + p cross- 
section scaled by the Glauber nuclear overlap function (Taa) 
in the centrality bin under consideration: 



Raa(pi 



cPNaa/ dprdy 
(Taa) ■ d 2 a pp /dp T dy ' 



(1) 



Deviations from Raa(pt) = 1 quantify the degree of depar- 
ture of the hard A+A yields from an incoherent superposition 
of AW collisions. Figure [2] compares the nuclear modifica- 
tion factor for r| in central (0-20%), semi-central (20-60%) 
and peripheral (60-92%) Au+Au reactions using the reference 
d 2 Opp/dprdy spectrum measured in p+p collisions iz7ll . As 
observed for high-/?;?- n° 1 8, 9], the r| yields are consistent with 
the expectation of independent AW scatterings in peripheral 
reactions (Raa ~ 1) but are increasingly reduced for smaller 
centralities. The pj dependence of Raa is flat above 4 GeV/c 
as seen also for the 71°. 

Figure [5] compares the Raa(pt) measured in Au+Au at 
= 200 GeV for r| (0-20% centrality), Jt° 0H and y Q 
(0-10% centralities). Whereas direct y are unsuppressed com- 
pared to the T^A-scaled reference given here by a next-to- 
leading-order (NLO) calculation |25, 30] that reproduces the 
PHENIX p+p photon data well Bill . %° and r| are suppressed 
by a similar factor of ~5 compared to the corresponding p+p 
cross-sections 1191 127T . Within the current uncertainties, light- 
quark mesons at RHIC show a flat suppression in the range pj 
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PHENIX Au+Au (central collisions): 
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FIG. 3: Raa{Pt) measured in central Au+Au at y/s NN = 200 GeV 
for n, 7t° 0-0] and direct y l25ll . The error bars include all point- 
to-point errors. The error bands at = 1 have the same meaning 
as in Fig. 13 T he baseline p+p — > y + X reference used is a NLO 
calculation 1 25, 30], that reproduces our own data well with 
theoretical uncertainties indicated by the dash-dotted lines around 
the points. The solid yellow curve is a parton energy loss prediction 
for a medium with density dN g jdy = 1 100 fl(h . 
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FIG. 4: n/ir ratio in Au+Au (centralities: 0-20%, 20-60%, 60-92%) 
compared to the ratio in p+p and d+Au at s /s^ = 200 GeV. The 
error bars include all point-to-point errors that do not cancel in the 
ratio of yields. The dashed curve is the PYTHIA 1 32] prediction for 
p+p at y/s = 200 GeV consistent with the asymptotic R^/rfs ~ 0.5 
measured in hadronic and e + e~ collisions in a wide range of cm. 
energies l27ll . 



= 4-14 GeV/c, independent of their mass (note that the r| is 
four times heavier than the ju°). The results are in agreement 
with expectations of in-medium non-Abelian energy loss of 
the parent parton prior to its fragmentation in the vacuum. The 
initial gluon densities needed to quench the high-pr hadrons 
by such an amount are of the order of dN g jdy = 1 100 (solid 
curve in Fig. [3} fill . 

An additional way to determine possible differences in the 
suppression pattern of Jt° and r| is to study the centrality de- 
pendence of the r|/ji° ratio in Au+Au collisions and com- 
pare it with the ratio in more elementary systems {e + e~ , p+p, 
d+Au). The r|/7i ratio in hadron-hadron, hadron-nucleus and 
nucleus-nucleus collisions is seen to increase rapidly with pj 
and flatten out above pj sw 2.5 GeV/c at an asymptotically 
constant R^/^o ~ 0.5 for all systems 1E7I1 . Likewise, in e + e~ 
at the Z pole (y/s = 91.2 GeV) one also finds R^/ n o ~ 0.5 for 

r| and K° at large scaled momenta x p = Phadmn/ Pbeam ^ 0.3 - 
0.7 12711 consistent with the range of fractional momenta (z) 
relevant for high-pr production discussed here. It is interest- 
ing to test if this ratio is modified in any way by final- and/or 
initial-state medium effects in Au+Au collisions at RHIC. 

Figure|4]shows R^^o {pr ) for three Au+Au centrality selec- 
tions and for p+p and d+Au collisions |27]. A fit to a constant 
for p T > 2 GeV/c gives R AuA f- 20 % = 0.40 ± 0.04, R dA «*f B = 

0.47 ± 0.03 and = 0.48 ± 0.03, where the quoted errors 
are the quadratic sum of statistical and systematic uncertain- 
ties. The Au+Au ratio is consistent within ~ la with both 
the essentially identical d+Au and p+p ratios. The R^/^o ra- 
tio shows thus no apparent collision system, centrality, o r p-r 
dependence. The dotted curve is the predicted PYTHIA 13211 
result for the p+p ratio at y/s = 200 GeV which is also co- 



incident with the world data measured in the same momen- 
tum range in hadronic, nuclear, and e + e~ collisions in a wide 
range of energies (y/s « 3 - 1800 GeV) JH. 

In summary, the transverse momentum spectra of r| mesons 
have been measured at mid-rapidity in the range pj = 2- 
10 GeV jc in Au+Au at ^JJ^ = 200 GeV. The invariant yields 
per nucleon-nucleon collision are increasingly depleted with 
centrality in comparison to p+p results at the same center-of- 
mass energy. The maximum suppression factor is ^5 in cen- 
tral Au+Au. The magnitude, pj, and centrality dependences 
of the suppression are the same for r| and 7t° suggesting that 
the production of light neutral mesons at large pj in nuclear 
collisions at RHIC is affected by the medium in the same way. 
The measured r|/7i ratio is flat with pj and amounts to R^/ K o 
= 0.40 ± 0.04. This value is consistent with the world value at 
high-^7- in hadronic and nuclear reactions and, at high x p , in 
e + e~ collisions. We conclude that all these observations are in 
agreement with a scenario where the parent parton first loses 
energy in the produced dense medium and then fragments into 
a leading meson in the vacuum according to the same proba- 
bilities that govern high-pr hadroproduction in more elemen- 
tary systems (p+p, e + e~). 
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